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We study a one-dimensional singular potential plus three types of regular interactions: constant electric field, harmonic oscillator 
and infinite square well. We use the Lippman-Schwinger Green function technique in order to search for the bound state energies. 
In the electric field case the unique bound state coincides with that found in an earlier study as the field is switched off. For non-zero 
field the ground state is shifted and positive energy "quasibound states" appear. For the harmonic oscillator we find a quantum phase 
transition of a novel type. This behavior does not occur in the corresponding case of an infinite square well and demonstrates the 
influence of quantum non-locality. 
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^ 1. Introduction 

i_i Singular potentials of Dirac 6 type are frequently used to de- 
scribe quantum systems, especially as toy models in many text- 

' pooks HI and in some papers on supersymmetry [2]. However, 
I there are a number of recent examples, based on 6 interactions, 

(f) of physical interest. To mention a few: Bose-Einstein conden- 

C*~) sation in a harmonic trap with a tight and deep "dimple" po- 
f ential, modeled by a Dirac 5 function is studied in . In J3l a 



singular part of the form 

W(x) = -a6(x) + bd'(x), a,b>0. 



(1) 



s£> 
O 



O 

> 



nonperturbative study of the entanglement of two directed poly- 
mers subjected to repulsive interactions given by a Dirac 5 func- 
Os tion potential has been carried out. Finally, in a series of three 
papers fl, light propagation in a one-dimensional realistic di- 
electric superlattice, modeled by means of a periodic array of 
k> Dirac delta functions, is investigated for the cases of transverse 
- electric, transverse magnetic, and omnidirectional polarization 
modes. Nevertheless, examples of the use of potentials related 
to the derivative of the Dirac delta function 5'{x) have appeared 
in the literature only very recently [6], and in a rather abstract 
context. Indeed, this may be due, in part, to the fact that there 
has been some controversy on the meaning of the 6'(x) potential 
|01, since different regularizations produce different reflection 
and transmission coefficients. 

In the present work we examine one-dimensional quantum 
Hamiltonians with potentials having a regular part Vq(x) and a 



Our efforts will concentrate on finding the bound states of these 
systems using the Green function as the fundamental tool. 

Let us consider a one-dimensional quantum system whose 
dynamical evolution can be governed by two different Hamilto- 
nian operators 

d 2 

H = --^ + V (JC), H = H + W(x), (2) 



dx 2 



the first one will be the regular free Hamiltonian Hq and the 
second one the perturbed Hamiltonian H (units have been cho- 
sen such that h = 2m = 1). Let us assume that Hq admits a 
Green function given by Gq(x, x', E) (the dimension of x and x' 
is here irrelevant). Then, i/r(x) represents a bound state of the 
total Hamiltonian H with energy E if and only if the Lippman- 
Schwinger bound state equation 



{x,x,E)W(x')tf/(x')dx . 



(3) 



'Corresponding author. 

Email addresses: larygaclarkson.edu (M.L. Glasser), 
manuelgadellal9gmail.com (M. Gadella), 
luismlSmetodos . f am. cie .uva. es (L.M. Nieto) 

URL: http : //metodos . f am. cie . uva. es/~ luismi/ (L.M. Nieto) 



is satisfied. In a previous publication (8|], this method was used 
to find the bound states in a one-dimensional system subject to 
a constant electric field plus a Dirac delta well. In this case, the 
free Hamiltonian was given by the kinetic term plus the elec- 
tric field, and the perturbation by the Dirac well. In the present 
paper, we shall extend this idea to the one-dimensional case in 
which the free Hamiltonian is a sum of the kinetic term plus 
either a constant electric field, a harmonic oscillator or an infi- 
nite square well. The perturbation is now given by the poten- 
tial W(x) introduced in Eq. ([TJ. To be explicit, we shall study 
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Hamiltonians of the form H - Hq + W(x) such that their unper- 
turbed part Hq has one of the following potentials: 



Voto = 
V (x) = 



—Fx , 

k 2 2 

16* 



F > 0, 



if 

oo if 



|jc| < c 
\x\ > c 



(4) 
(5) 

(6) 



We shall see that the presence of the 5'{x) term has a dramatic 
effect in the composition of the bound states of Hq. In all cases, 
the presence of the perturbation term W(x) in the potential, may 
create difficulties in interpretation. These can be overcome by 
the specification of matching conditions for wave functions at 
the origin, which is the singular point. This is the approach 
followed by Kurasov and coworkers J^lOU; it avoids the con- 
troversy on the meaning of the 6'(x) potential [ 7] and gives rea- 
sonable physical results |6]. These matching conditions show 
that the wave functions and their derivatives may be discontin- 
uous at the origin. This is opposed to the usual procedure when 
solving quantum mechanical problems, because the continuity 
of the wave function is normally taken for granted, a fact that 
must be discarded here. 

The presence of Six) and 6'{x) together in the potential im- 
plies that the Schrodinger equation necessarily be a distribu- 
tional equation. In addition as the wave functions are discon- 
tinuous at the origin, we have to define the product of a discon- 
tinuous function if/(x) times Six) and 6'(x). This can be done as 
follows 101: 



u w ^ «A(0+) + *A(0-) s , . 
\p{x)o{x) = 5(x) . 



(7) 



tf/(x)6 (x) = 6 (x) 6(x) , 



where <p(0+) and <p(0-) are the right and left limits of <p(x) at 
the origin (here, <p{x) is either tfr(x) or its derivative). The two 
equations © are essential in the present study. 

The next section consists of three parts. First of all, starting 
from the explicit form of d3j with W(x) - -aS(x) + bd'(x), 
we obtain a determinant whose roots give the energy values 
for the bound state of the total Hamiltonian H — Hq + W(x). 
Then, in the following three subsections we search for bound 
states for the total Hamiltonian, where Hq is given by (0]i (con- 
stant electric field), (0 (harmonic oscillator) and © (infinite 
square well) respectively. The most striking result appears in 
the case (0). For values of the parameter b 2 above a threshold, 
the harmonic oscillator loses its bound states which become res- 
onances grouped in complex conjugate pairs. If b 1 takes values 
below the threshold, there is an infinite number of bound states. 
In the case ©, there exists the possibility for a unique bound 
state with negative energy in a narrow range of values of b. This 
negative value has no lower bound. The results are summarized 
in the Concluding Remarks. 



2. The search for bound states 

We can apply formula © to all the three cases considered 
above, because for them we explicitly know the Green function 
Gq(x,x',E) of the corresponding Hq. In this case, (01 has the 
form 



= Jg (x,x', 



E)[-a6(x') + bS'(x')] iff(x') dx , (8) 



In all other cases, we shall use the Green function G f ree (x, x', E) 
for the Hamiltonian of a free particle, and then apply Q with 
W(x) as the total potential. 

From now on we shall use the following notation for simplic- 
ity 



Q'+j 

dx-dx'j 



Gq(x,x',E), i,j = 0, 1. 



(9) 



Then, from <£3j we readily obtain 

2(A(x) = -aG (x,0)[(A(0+) + (/r(0-)] 

-b[Gl\x, 0)(/<0+) + G (x, 0)i/r'(0+) 



+G n (x, 0)i/<0-) + G(x, 0)^(0-)] . (10) 



Next, we introduce the following notation to simplify our ex- 
pressions: 



A := G (0+,0) = Go(0-,0), 
B := GJ^O-kO) = G^°(0-,0), 
C := G^(0-,0) = G l \0+,0), 
D := Gi'(0+,0) = G ( V(0-,0), 



(11) 
(12) 

(13) 
(14) 



where the signs plus and minus after a denote right and left 
limits, respectively. Identities such as Go(0+,0) = Go(0-,0) 
are intrinsic properties of the Green functions considered in this 
paper. 

If we take right and left limits at the origin in fllOl . we obtain 
respectively, 



a b 
1+-A+-B 

2 2 



<A(0+) + 



a b 
2 2 



<A(0~) 



+ ^A.A'(0+) + ^A ( A'(0-) = 0, (15) 



a . b „ 




2 2 


<A(0+) + 



a . b „ 
1+-A+-C 
2 2 



<A(0-) 



+^A l A'(0+) + ^A ( A'(0-) = 0. (16) 

When we perform the same operations with the derivative of 
( [Tol l, we find 



a „ b 
-C+-D 

2 2 



a b 
-B+-D 

2 2 



<A(0+) + 



<A(0+) + 



a b 
-C+-D 

2 2 



<K0-) 



iA'(0+)+-C<A'(0-) = 0, (17) 



a b 
-B+-D 

2 2 



b 

1 + -B 

2 



•K0-) 

«A'(0-) = 0. (18) 



2 



Taken together, the four linear equations (Tl~5T>— (fT8l> give a 4 x 4 
system which has a non-trivial solution if and only if 



det 



This determinant can be simplified by subtracting the third col- 
umn in from the last, adding the fourth row to the third and then 
subtracting the second row from the first. Then, by subtracting 
the second column from the first one obtains, 



1 + § A + f B 


§ A + | B 


§A 


-A ^ 

2 71 




1 + § A + §C 


§A 


|A 


- C + - D 

2 2 


|C+|D 


1 + |C 


|C 


§B+§Z) 


§B+§D 


§* 


1 + §B , 



det 



-l-f(B-C) 
1 + f A + |C 







§(B + C) + M) 1 + |(B + C) 



(19) 



Equation (fT9l will give us the values of the energy for the bound 
states in three cases under consideration, which we analyze sep- 
arately below. 

2.1. Constant electric field 

In this case, Vo is given by |@). The appropriate Green func- 
tion Gq(x, x' , E) := Go(x, x'), given in (8J, is 



Go(x,x) = -—Ai\ p _ U3 



Since £ is fixed, the notation Gq(x, x') will be used where no 
confusion is possible (see (0). In (EOt . Ai(x) and Bi(x) are the 
Airy functions of first and second kind, respectively, and x > and 
x< are the maximum and minimum of x and x' respectively. As 
we already pointed out in (fTTT i. this Green function is contin- 
uous at the origin with respect both arguments, so that, if we 
denote z = -E/F 2 ^, 



A = G (0+,0) = Go(0-,0) 



FV3 



Ai(z)Bi(z). 



(21) 



At the origin, partial derivatives with respect to the arguments 
x and x' have jumps and fulfil the identities given in (fT2l> — (TT~4b. 
Then, we have 



B = -7TAi(z)Bi'( Z ), 

C = -7rAi'(z)Bi(z), 

D = -nF lji Ai'(z)Bi'(z). 



(22) 
(23) 
(24) 



By inserting these values into ( fT9l ), we obtain a transcendental 
equation that gives the energy values for the total Hamiltonian 
©// = //o + W{x), with V (x) as in ©. A plot of the bound 
states resulting from the transcendental equation ( fT9l . which in- 
volves Airy functions, is shown in Figure 1 for different values 
of a and b. 

If we take b — in the perturbation part of the potential W(x) 
given by ([TJ, then Eq. (fT9l ) reduces to 1 + aA = 0, which is 
the problem treated in On taking the limit F — > we find 



6 . 

4 

2 

Kl - 

-2 - 
-4 - 
-6 - 
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Figure 1: Bound states for the electric field F = 1, resulting from 
Eq. ([19} with $2l\-$24[ and for b = -1 (dotted curve), Z? = (solid 
curve), and b = 1 (dashed curve). Units are such that h = 2m = 1. 



that the energy is E = -a 2 /4 as expected. In the case of low 
field strength, where F is small, the properties of Airy functions 
yield, 



A = — 



1 



D = 



(25) 



For 5 s 0, (O simplifies to aA = b 2 AD - 1. Using (|25]) and 
taking into account the definition of z, we obtain the value of 
the energy of the bound state in the limit F — > 0, which is 



E = - 



4a 2 



(4 + fc 2 ) 2 ' 



(26) 



in agreement with [6] (remember that throughout the present 
work we have chosen m =1/2). 

A striking result found in [8] is the existence of positive en- 
ergy solutions to the bound state equation for the 5 well, which 
correspond to "quasibound" levels embedded in the continuum 
of scattering states. To check whether these are seriously af- 
fected by the presence of the 5' potential we examined the evo- 
lution of the spectrum, numerically for the case a = b = 1 and 
m = 1/2, as F increases. The findings, shown in Figure 2, sug- 
gest that the quasibound state spectrum is similar to that found 
in ||8J]. In particular, the ionization field F c = 0.076 I8[] is de- 
creased by about 10% to 0.062 due to the influence of b + 0. 

2.2. The Harmonic Oscillator 

The next Hamiltonian to be considered is the one coming 
from the potential in Eq. (01: 



H 



dx 2 



ad(x) + b6'(x) 



k 2 x 2 



(27) 



where the total potential is a sum of a harmonic oscillator well, 
where the form of the coupling constant has been chosen for 



3 



0.02 0.04 0.06 



0.10 - 



0.05 J 



ftj 0.00 



-0.05 - 



-0.10 - 




Note that d32l is strictly positive for k > and E < (since 
> whenever x > 0). This shows that the action of the 
potential given in <(5j removes the negative energy bound state 
that appeared with the singular potential and a > 0. 



4 - 



2 - 



0.02 0.04 0.06 



Figure 2: (Color online) Bound states of the constant electric field case 
for a = b = 1 (units are h = 2m = 1). Solutions appear only in a 
small interval of the electric field: < F < 0.06151359. Note the 
presence of an infinite number of loops in that region, which in the 
figure appear distorted due to the numerical approximation. In the 
inset, an amplification of this area is shown. 



convenience, and the singular potential at the origin. By setting 
k — 0, i.e., dropping the oscillator term in the potential, we can 
verify that our method gives the same value as obtained in 01 . 
It is well known that the Green function for the free particle is 
given by (2m = K = 1) 



Go(x, x', E) 
If we use 
A = 



exp [iy[E\x - x'\] 



E > 0. 



2/ VI 

in (HU)-® with a + + b, k = 0, we obtain 



(28) 



1 



D = — V#. 

2 



2iVI' 

Then, by using ( 1291 in (fT9l l. we get the equation: 
1 



1 + 



b 2 

2; VI 4 



(29) 



(30) 



Eq. (l30l makes sense if and only if E takes negative values. 
Indeed, this is what happens and we find again the result given 
in Eq. d26l l. which coincides with the one given in 

Now let us assume that the three coefficients in (l27l i are dif- 
ferent from zero. In this case, one has to calculate the Green 
function corresponding to Vq in (0 which has been obtained in 
| n|] and gives 



A = 



rft-f) 



2V£r(i-f) 

Accordingly, (T% becomes: 

b ^(l-l) a r(i-f) 
2 r 2 (|-f) + VSr(|-f) 



B = C = , D = -M 



+ 2 = 0. 



(31) 



(32) 



-2 



-1 



Figure 3: Plot of the bound states from Eq. d32t as a function of b (units 
are k = h = 2m = 1), for a = 1 (dotted curve), a = 2 (solid curve), and 
a = 4 (dashed curve). 



Concerning the oscillator bound states, note that the first term 
on the left hand side of ( l32l is always positive. The second term 
could be either positive or negative, but for large values of k one 
expects that eventually this term will be smaller than 2. In this 
case, the left hand side of d32l ) is positive and, therefore, no 
solutions for E exist; i.e. no bound states exist for d27l >. To 
show when this happens, let us write ( f32b as 



u2 2 

b a 



+ 2 a + 4 



, where a : - 



r (|-f) 

r(i-f) 



(33) 



The discriminant of (1331 ) is ^(a 2 /k) - 4b 2 ; if it is negative, 
equation d33l has no real root a and therefore, no solutions E 
exist when 



b 2 > a 2 /{4k) . 



(34) 



Hence, there are no solutions E for Eq. (l32l for values of 
b above a/(2*\[k) (or below —a/ (2 Vk), since this equation is 
symmetric with respect to a change of sign in b), and thus there 
are no bound states! This is a very striking result, that can 
be appreciated in Figure 3: the action of the singular potential 
W(x) = -a5{x) + b6'(x) erases all bound states of the harmonic 
oscillator well provided that the absolute value of b is bigger 
than a critical value. This effect is produced solely by the pres- 
ence of the 6'(x) interaction with sufficient intensity. For val- 
ues of b in the interval (-a/(2 yfk), a/(2 yfk)), including b — 0, 
equation OTb has an infinite number of solutions grouped in 
pairs (due to the double sign in the solution a in (l33l l. see (|35T>). 
This is illustrated in Figures 3 and 4. 




- 



2 4 6 8 10 

a 

Figure 4: Plot of the bound states from Eq. i32l as a function of a, (units 
are k = h = 2m = 1), for b = (dotted curve), b = 1 (solid curve), and 
b = 3 (dashed curve). 



that for each value of b ( f37l > has an infinite number of positive 
solutions, each giving the energy of a bound state. However, 
for values of b in the interval [2 y/n, a/2 + 4-n] there exists an 
additional negative solution for E. This energy is zero for b = 
a/2 + 4-k, but goes to -oo as b approaches to 2 \fn from above. 
In Figure 6 we display the bound state spectrum as a function 
of a for b = 1,5,10. The vertical line E — n 2 comes from the 
numerical algorithm and is spurious. 



b 




E 



-50 50 



For values of b in the range given by d34} equation d32l has 
complex solutions. Solving d33l for a, we obtain 



a = 



r(I-f) 



l 

b 2 



a 2 - 

- -4b 2 
k 



(35) 



This gives for b, satisfying (I34} . complex solutions of (l32l for 
£ grouped in pairs. Since the gamma function has the property 
Hz*) = [r(z)]*, where the star denotes complex conjugation, 
these pairs are complex conjugates, so that if E^+iy is a solution 
of d35l >. then Er - iy is also a solution. These solutions can be 
considered to be resonances. Unlike the case b — 0, treated in 
Jst], both the even and odd ascillator states are affected by the 
perturbation. 

2.3. The infinite square well 

The last case under our study is the Hamiltonian H = Hq + 
W(x), where Vo is given by (O. The corresponding Green func- 
tion has already been calculated in 01211 and gives for A,B,C 
and D in Eqs. (fTTT> — <TT~4-b the following values 



tan(c VE) 

27T VE 



B = C = 0, 



D 



2 tan(c VF) 



(36) 



2c being the well width of the potential ©. If we insert these 
values into (fT9] >. we obtain the transcendental equation 



tan(c VE) b 2 - 4n 



2a 



which has an infinite number of solutions. Note that if a — 
this equation is meaningless. 

In order to search for solutions of (|37| |. let us choose c = 1, 
fix some value of a, and then plot the resulting function of b in 
terms of the energy. This is done in Figure 5, and we can see 



Figure 5: (Color online) Plot of the bound states from Eq. {37} as a 
function of b (units are c = h = 2m = 1), for a = 1 (dotted curve), 
a = 5 (solid curve), and a = 10 (dashed curve). 



2 - 



10 



15 




Figure 6: (Color online) Plot of the bound states from Eq. {37} as a 
function of a, (units are c = h = 2m = 1), for b = 1 (dotted curve), 
b = 5 (solid curve), and b = 10 (dashed curve). The vertical line 
at E = n 2 does not correspond to any bound state: it is a spurious 
solution coming from the numerical calculation. Observe that there is 
an infinite number of bound states (with similar structure) for E > 0, 
but not for E < 0. 



(37) 3. Concluding remarks 



To summarize this work, the Lippman-Schwinger Green 
function technique has been used to study the energy of pos- 
sible bound states for a quantum Hamiltonian of the form 
H = --fj-? + Vq(x) - ad(x) + bd'(x). Three particular cases have 
been studied. 



5 



(i) Vq(x) = -Fx with F constant. The case b — was studied 
previously. In the limit F — > 0, we obtain one negative eigen- 
value which coincides with the value found in [6] that uses a 
quite different technique. 

(ii) Vq(x) = k 2 x 2 /l6. In this case, there can be no state with 
negative energy. But the most striking result is the existence of 
a quantum phase transition El characterized by a threshold for 
the absolute value of b (that also depends on a and k). Below 
this threshold, there is an infinite number of bound states which 
group into pairs. Above this threshold, the infinite ladder of 
harmonic oscillator states evaporates and is replaced by an infi- 
nite number of resonances, which group in complex conjugate 
pairs. 

(iii) If Vq(x) is an infinite square well potential, the singular 
interaction W(x) does not eliminate the bound states of the free 
Hamiltonian Hq = -cfi/dx 2 + Vo(x), which are merely shifted. 
For values of b in a narrow interval, there exists a unique nega- 
tive energy bound state whose absolute value can be arbitrarily 
large. 

The difference in behavior between the square well and the 
harmonic well cases can only be due to the influence of bound- 
ary conditions, "soft" in the harmonic case and impenetrable in 
the square well, in a region where the singular potential is zero. 
This provides a striking demonstration of non-locality. 

Some other examples of more realistic potentials Vq(x) 
should be studied to obtain more conclusions about the emer- 
gence of the quantum phase transitions mentioned above. To 
analyze this phenomenon in two and three dimensions is also a 
very attractive project, that is presently in progress. 
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